Phase transition without global ordering in a hierarchical scale-free network 
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(N ■ 

I ' We study the site-bond percolation on a hierarchical scale-free network, namely the decorated 

, (2,2)-flower, by using the renormalization group technique. The phase diagram essentially depends 

■ on the fraction of occupied sites. Surprisingly, when each site is unoccupied even with a small 

probability, the system permits neither the percolating phase nor the nonpercolating pheise but 
rather only critical phases. Although the order parameter always remains zero, a transition still 
exists between the critical phases that is characterized by the value of the fractal exponent, which 
measures the degree of criticality; the system changes from one critical state to another with the 
jump of the fractal exponent at the transition point. The phase boundary depends on the fraction 
of occupied sites. When the fraction of unoccupied sites exceeds a certain value, the transition line 
between the critical phases disappears, and a unique critical phase remains. 

^ PACS numbers: 89.75.Hc 64.60.aq 89.65.-s 

Hierarchical lattices have occupied a prominent position in the field of statistical physics, because the cooperative 
behaviors on such lattices are exactly solved by the renormalization group (RG) technique. Recently, hierarchical 
' lattices have also been applied in the context of complex networks [H-lEl- Adjusting the construction rule enables us to 
^ ' generate exactly solvable networks having properties common to real networks: a scale-free (SF) degree distribution 
I ■ P{k) oc {k denotes degree), a mean shortest path length with logarithmic dependence I ^ 0(ln A^) (A^ being the 
' ^ ■ number of sites) , and a highly clustering coefficient C @ . 

Dynamics (such as percolations and spin systems) on hierarchical SF networks provide an impetus to further 
investigate the relation between network topology and dynamics because of the abnormal phase transitions on such 
networks. Hinczewski and Berker Q analyzed the Ising model on a hierarchical SF network, called the decorated 
(2,2)-flower (Fig.[T]), to show that the system undergoes an inverted Berezinskii-Kosterlitz-Thouless (BKT) transition, 
] which means that the system shows a BKT singularity above the transition temperature. A similar abnormal behavior 
' holds for the case of bond percolation 0, |^. The bond percolation on the decorated (2,2)-flower has the percolating 
] phase, which contains a unique giant component whose size is of the order 0(N), and the critical phase (also known 
QO . as the partially ordered phase or the "patchy" phase for the Hanoi network ^^). In terms of RG, the critical 
' phase is characterized by the RG flow converging onto the line of the nontrivial stable fixed point in contrast to the 
; nonpercolating phase containing only finite size clusters. 
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FIG. 1: (Color online) Realization of the decorated (2,2)-fiower F„ (and _F„) with n = 0, 1, 2, 3. The root sites are indicated 
by the open circles. Fn (solid black lines) is recursively constructed (see the text), and Fn is obtained by adding the shortcuts 
(red dashed lines) to Fn. Note that the shortcuts remain the same in each iteration. 
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FIG. 2: (Color online) Phase diagrams for several values of p,,: (a) ps = 1.0, (b) ps = 0.9, (c) ps = 53/64, and (d) Pa = 3/4. 
The bold black line indicates the phase boundary. The blue line at p = 1 indicates the percolating phase, and the red line at 
p = indicates the nonpercolating phase. Each inset shows the lines of stable and unstable fixed points. The solid blue and 
red dashed lines indicate p*{p,ps) and p**{j>,ps) (stable fixed point) and pc{p,Ps) (unstable fixed point), respectively. 



In the critical phase, where the order parameter is zero, the system is characterized by the fractal exponent tp 
The fractal exponent is defined as the logarithmic derivative of the mean size S'inax(^) of the largest cluster, 
ijj = limAr_^oo log^r >S'inax(-^), which mimics df /d for d-dimensional Euclidean lattice systems, where df is the fractal 
dimension of the largest clusters. In the critical phase tp takes a positive value less than one, < tp < 1, while ip = 1 
indicates the percolating phase and tp = the nonpercolating phase. The mean number rig of clusters with size s per 
site (or the cluster size distribution, in short) obeys the power-law Us cx s~'^ with varying exponent r = 1 + tp~^ over 
the entire region of the critical phase. Such a critical phase is also observed in growing random networks (see and 
references therein). 

Previous studies revealed that the phase diagram of dynamics on complex networks is not simply described by 
globally ordered state and disordered state. The concept of partially ordering may be needed for the complete 
understanding of the dynamics on the complex networks. In this paper, we investigate the site-bond percolation on 
the decorated (2,2)-flower to demonstrate the existence of a hidden phase transition without globally ordered state. 
In this network, defects of sites may have significantly different effects from those of the defects of bonds, because 
the defect of a site with a high degree is directly connected to the removal of infinitely many bonds, and also, the 
sites with high degrees are interconnected. We demonstrate that the site dilution dramatically changes the phase 
diagram. Surprisingly, even if the fraction of the removed sites is very small, global ordering is broken, and only 
the partially ordered state is permitted. Furthermore, although the order parameter is zero in the entire region, a 
transition between two different critical phases still exists, which is characterized by a discontinuous jump of the 
fractal exponent. When the fraction of the removed sites exceeds a certain value, the boundary between the two 
critical phases disappears, and a unique critical phase remains. 

A special class of hierarchical SF networks, the (decorated) (m, 'i;)-flower, was introduced in The (2,2)-flower Fn 
of the n-th generation is recursively constructed as follows (Fig. [IJ: At 71 = 0, Fq consists of two sites connected by 
a bond. We call these sites root sites. For n > 1, Fn is obtained from Fn^i, such that each existing bond in Fn-i is 
replaced by two parallel paths consisting of two bonds each. As illustrated in Fig. [U the decorated (2,2)-flower Fn 
of the nth generation is given by adding the shortcuts (red dashed lines) to Fn- The number of sites iV„ of Fn is 
Nn = 2(4" + 2)/3, and F„ has P{k) oc k'^, ioclnN, and C - 0.820 0. 

Let us consider the site-bond percolation on the decorated (2,2)-flower. Each site is occupied with the probability 
Ps. The open-bond probability of bonds constructing Fn is p, and that of the shortcuts is p. The phase diagram of the 
model is obtained by the RG technique 0,^01. Let P^"^ be the probabihty that both roots are in the same cluster 
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FIG. 3: (Color online) (a) Flow diagram for ps = 0.9. Arrows indicate the flow of Eq. IT}. The solid blue and red dashed lines 
indicate the stable and unstable fixed points, respectively. The inset shows the fractal exponent as a. function of p. Here 
■i/) = 0.6504, 0.9457, and 0.9994 for p < pc(p,Ps),p = Pc{p,Ps), and p > pc(p,Pa), respectively, (b) AP'"' as a function of P'-"K 
The open and filled circles indicate the unstable and stable fixed points, respectively. In the inset of (a) and also in (b), we set 
Pa = 0.9 and p = 0.1 [the horizontal dotted line in the main panel of (a)]. 



after a site-bond percolation trial on Fn- Here we do not inquire whether the two roots are occupied or not. In the 
large size limit, the system is regarded as percolating if the percolation probability P = lim„^oo P'"^ is nonzero 
We find the recursive equation of P'"-' as 



p{n+l) 



= P 



2p,(P("))2-p2(p(n))4l ^ (1) 



where q = 1 — p, and the initial value is set to P^^^ = p. The fixed point P = P^") = p("+i) of Eq. ([IJ satisfies 

piP) = l- ,^ P2,2 - (2) 

(1 -psP^y 

We obtain the phase diagram from the RG technique [Eq. ([T|)] and the fixed-point condition [Eq^^©]. 

The phase diagram for ps = 1, i.e., for the bond percolation with p and p, is already known 0, Q [Fig. [2ja)]. For 
a fixed p [Q < p < pc = 5/32). there are two stable fixed points: < P = p*{p,ps) < 1 and P = 1, which correspond 
to the critical and percolating phases, respectively. Furthermore, one unstable fixed point between these two fixed 
points gives the phase boundary: P = pc{p,Ps)- For p > pc, there is only one stable fixed point at P = 1, so that the 
system is always percolating. 

When a site dilution is added (3/4 < ps < 1), the phase diagram dramatically changes [Fig. d^b)]. The flow 
diagram with ps = 0.9 is shown in Fig. ^a). For a given p {0 < p < pc), there are two nontrivial stable fixed points at 
P = p*{p,Ps),p**{p,Ps) and one unstable fixed point at P = pc{p,Ps), where < p*{p,ps) < Pc{p,Ps) < P**ip,Ps) < 1- 
Here the (in) stabilities of the fixed points are given by AP(") ee p("+i) - p(") [Fig. gjjb)]. We find that all RG flows 
converge onto nontrivial flxcd points (0 < P < 1). The system has neither the percolating phase nor the nonpercolating 
phase, but rather, it has two different critical phases, which we call critical phases I and II. In the entire region, Us 
obeys a power-law with a corresponding exponent (not shown). Also, the mean cluster size (s) = '^gs'^ris always 
diverges, because ip > 1/2, which means r < 3 (see below). For p > pc, one stable flxed point P = p**{p,ps) exists, 
and the system is always in the same critical phase. 

A kind of phase transition between critical phases may exist when we pass an unstable fixed point, even though 
the order parameter remains zero in the entire parameter region. Let us calculate the fractal exponent of the mean 
size S'root(A^) of the cluster to which a root belongs, defined a.s ip = lim„^oo dlnS'ioot(A^n)/dln A'n. We consider 
three quantities on P„: (1) the probability t^^\p) that both roots are connected to the same cluster of size fc, (2) 
the probability s\l^i (j>) that the left (right) root is connected to a cluster of size k (l) but that these clusters are 
not the same, and (3) the mean number u^"-* (p) of clusters of size k to which neither of the roots is connected. For 
the sake of convenience, the roots are not counted in the cluster size k or I for t[,"''(p) and s^"^(p). We introduce 
the generating functions r„(a;), Sn{x,y), and Un{x) for ^kjip)- ^^^d u[."^(p), where T„(a;) = I]fclo 

Sn (x, y) = J^kLo E/^o ^k^i (p)^'^y' J ^i^d Un (x) = J2T=o "fe"^ ip)x'' . The self-similar structure of P„ allows us to obtain 
recursion relations for these generating functions: 
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FIG. 4: Fractal exponent -0 as a function of p for several values of ps'. (a) ps — 1.0, (b) ps — 0.9, (c) ps — 53/64, and (d) 
Ps = 3/4. 



Tn+iix) ^ {pIx^ + 2psqsx)T^{x) + Aplx'^Tl{x)Sn[x,x) + 4psqsxT^{x)Snil, x) + 2psxT^{x)S„{x, l)S'„(l,x), (3) 
Sn+i{x,y) = {S'„(x, l)S'„(l,y) +Ps'S'„(x,y) [a;T„(x) + yT„(?/)] + [r„(a;)T„(2/) + T„(a;)5„(l, y) + T„(t/)S'„(a;, l)]}f4,) 
C/„+i(x) = 4[/„(x)+2{p,x5„(l,a:)5„(a:,l) + (7, [5„(l,.T) + 5„(a;,l)]}, (5) 

where qs = 1 — Ps- The corresponding generating functions Tn{x), Snix,y), and Un{x) on F„ are given by T„(.t) = 
T„{x) + pSnix, x), Sn{x,y) = qSnix,y), and Un(x) = C/„(x), respectively. Here we denote the mean fractions of 
the clusters including two roots and either of the roots by t„ = ^Tn{x)\x=i and cf„ 
respectively. Similarly to @, we obtain the recursive equation for t„ and cr„ as 



7V-iA5„(x,x)U=i, 



fri+l 
Tn+l 



_( 2g(l-p,(p("))2)(l+p,p(")) 4g(l-p,(P("))2)(l-p,p(»)) ^ 

_i \2{l +p,P^^)) [1 - g(l -p,(p("))2)] 4 [1 - g(l -p,(p("))2)(i -p^pH)]^ 

(«)g(n)(l _p^(p("))2) \ 



"iV„+i ^2p,p(") [2-P(") -29Q(")(l-p,(p("))2)] 

1 Q(i+P.P) Q(i-P.P) ^ 

1 P{l+p,P){l-p,P) _ Q(l-p.P) 
^2 l-PsP'^ l-p,P^ 




for n 2> 1, 



(6) 
(7) 



where Q = 1 — P and Q^"' = 1 — P*^"', and we used the fixed-point condition [Eq. ([2])]. Then, the largest eigenvalue 
A of Eq. ([7]) gives the fractal exponent as V' = 1 + In A/ In 4. The inset of Fig. [3] (a) shows the fractal exponent -0 for 
Ps = 0.9 and p = 0.1. When p increases, the fractal exponent jumps from a nonzero value ip = 0.6504 to another 
larger value ip = 0.9994 at the transition point pc{p,Ps)- 

The phase diagram and the corresponding V-'-curvc depends on ps (Figs. [2] and S]). When ps decreases to below 
Ps — 3/4, the unstable fixed point disappears, and the two lines of the stable fixed points merge. Then ^ increases 
monotonically with p, but never jumps (Fig. [4]). For ps < 3/4, the boundary between critical phases I and II 
disappears, and the unique critical phase remains. 

Finally, we mention the phase diagram for p = 1 (Fig. [5]). Figure [S] also indicates that the percolating phase is 
destroyed by infinitesimal site dilution. This destruction is related to the fact that the sites with high degrees are 
interconnected in a hierarchical manner. Moreover, the phase boundary between critical phases I and II disappears 
at {p,Ps) = (1/4, 3/4). This resembles the phase diagram of a vapor-liquid transition. 

To summarize, we have investigated the site-bond percolation on the decorated (2,2)-flower. The phase diagram 
essentially depends on the fraction of the occupied sites ps'. the system has the percolating and critical phases at 
Ps ~ 1, two critical phases (I, II) for 3/4 < p,, < 1, and a unique critical phase for ps < 3/4. We suggest that 
critical phase I is caused by shortcut insertion, while critical phase II is due to the site dilution. The shortcuts of 
the decorated (2,2)-flower, which connect high-degree nodes, easily give rise to a partially ordered state. Further, site 
dilutions in such connections immediately break a globally ordered state into a partially ordered state. 

We demonstrated the existence of a new type of phase transition without global ordering, namely, a transition 
between two critical phases. Our result seems to hold for the other decorated {u, w)-flower (u, v > 1), although it is an 
open question whether such a transition exists on other complex network models. Site dilutions in complex networks 
may also have extreme effects on other dynamics (e.g., spin systems and the contact process). We hope that further 
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FIG. 5: (Color online) Phase diagram of the case with p = 1. The bold black line indicates the transition line between the 
critical phases. The blue line (at ps = 1) and the red line (at p = 0) indicate the percolating and nonpercolating phases, 
respectively. 



studies on phase transition without global ordering give a new perspective on the critical phenomena emerging in 
complex networks. 
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